Semiflexible Polymer in a Strip 
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We study the thermodynamic properties of a semiflexible polymer confined inside strips of widths 
L < 9 defined on a square lattice. The polymer is modeled as a self-avoiding walk and a short-range 
interaction between the monomers and the walls is included through an energy e associated to each 
monomer placed on one of the walls. Also, an energy et, is associated to each elementary bend of the 
walk. The free energy of the model is obtained exactly through a transfer matrix formalism. The 
profile of monomer density and the force on the walls are obtained. We notice that as eb is decreased, 
the range of values for e for which the density profile is neither convex nor concave increases, and 
for sufficiently attracting walls (e < 0) we find that in general the attractive force is maximum for 
£6 < 0, that is, for situations where the bends are favored. 
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I. INTRODUCTION 

Polymers are often modeled as self- and mutually 
avoiding walks placed on a lattice, and much has been 
learned about their thermodynamic properties through 
such models HJ^]. Grand-canonical models of this kind, 
where the number of monomers incorporated into the 
chain is allowed to fluctuate controlled by an activity 
z = cxp(/i/fcsT), display a phase transition at some 
value of the activity (for infinite chains, that is, in the 
polymer limit). This transition is discontinuous in the 
one dimensional case d = 1 ||, and continuous for d > 2. 
Rather precise estimates of the critical value of z were 
found in two dimensions through transfer matrix calcu- 
lations H and series expansions ||. Also, exact values 
for the critical exponents are available in this case || . 

More recently, properties of such models on lattices 
limited by walls have attracted much interest 0, fol- 
lowing studies of magnetic models in the same situation 
||. The short range interaction between the wall and 
the polymer may be introduced associating a Boltzmann 
factor lu = exp(— e/ksT) to each monomer placed on 
the wall, so that e < corresponds to attracting walls 
while repulsive walls are described by e > 0. The grand- 
canonical partition function for a model of a single chain 
is 



Y(x.u) 



(1) 



where N is the number of monomers in the chain, the 
sum is over all configurations of the chain with the initial 
monomer placed on the wall, and N w < N stands for the 
number of monomers located on the wall. Such a model 
shows interesting features, and even in the limit where 
the self avoidance constraint is neglected (the so-called 
ideal chains) one finds that for sufficiently large values of 



lu > loq the surface polymerization transition will occur 
at a lower value of the activity z than the one in the 
bulk g|. The point (zo,coo) m the phase diagram where 
the bulk (also called ordinary) transition line meets the 
surface transition line is called the adsorption transition 
point. In two dimensions, such models have been stud- 
ied through transfer matrix calculations |l0| and series 
expansions jfj]]. Additional walls may be added, confin- 
ing the polymer inside a strip, slab or pore Q, and the 
force applied on the walls in such situations is of interest 
even from the point of view of applications of polymers 
as adhesives M] . The model of ideal chains confined in a 
slab has been studied in the past [|l2|,[l3) , and it was found 
that the force on the walls is attractive if u> exceeds the 
adsorption value uiq. In the case of self-avoiding chains 
confined in a strip on the square lattice, transfer ma- 
trix calculations show that attractive forces appear for 
lu below luq 1 13 1 . Also, the profile of the monomer den- 
sity inside the strip was obtained G3], and, unlike what 
happens for ideal chains, for a self- avoiding chain pro- 
files which are neither convex nor concave are found for 
a range of values of lu. Finally, similar techniques have 
also been applied to shed light on the scaling behavior of 
such models for strips with periodic boundary conditions 
(cylinders) as the width of the strips becomes large [jl5| . 

Another generalization of the original polymer model 
is to introduce an energy associated to the operation of 
bending the chain. On hypercubic lattices, the elemen- 
tary bends will always be at right angles, and an energy 
eb may be associated to each of them. This semiflexible 
polymer problem (also called persistent or biased walks) , 
has been studied some time ago 16-18 . Recently, this 



model has attracted renewed interest , since it may de- 
scribe some relevant aspects in the protein folding prob- 
lem [^9|. The thermodynamic properties of the model 
have been studied on the Bethe lattice 20 and, equiva- 



1 



lently, in the Bethe approximation |2l]]. The end-to-end 
distance of semi-flexible chains on Bethe and Husimi lat- 
tices was obtained 22 1. In this paper we study the ther- 



modynamic behavior of a semiflexible polymer confined 
inside a strip. The partition function of the model may 
be written as 



z to u>i , 



(2) 



where Nf, is the number of elementary bends in the con- 
figuration and LOb — exp(— e^/fe^T) is the Boltzmann fac- 
tor associated to each elementary bend and the sum is 
over all configurations of the chain. We define a transfer 
matrix for the model and obtain the grand-canonical par- 
tition function in the thermodynamic limit, determined 
by the largest eigenvalue of this matrix. In order to be 
able to obtain the distribution of monomers in the trans- 
verse section of the strip, we define a position dependent 
activity for the monomers. Also, the force applied by the 
polymer on the walls is calculated, as a function of the 
width of the strip, u> and Wf,. All thermodynamic proper- 
ties are calculated at the polymerization transition, that 
is, at the value of the activity of a monomer for which 
the number of monomers incorporated into the polymer 
diverges. 

In section [jl] the model is defined in detail and its so- 
lution is presented. The results we find for the thermo- 
dynamic behavior of the model are shown in section III 
as well as final conclusions and discussions. 



II. DEFINITION OF THE MODEL AND ITS 
SOLUTION 



number of activities n a (m) is equal to m/2 + 1 for even 
values of m and (m + l)/2 for odd values of m. 
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FIG. 1. a) The five connectivity configurations for m = 2. 
b) Portion of the chain, with the number of each connectivity 
configuration indicated 



The self avoiding chain is constrained inside a strip of 
width m defined on the square lattice in the (x, y) plane, 
so that < x < m. The chain runs through the whole 
strip, from y —> — oo to y — ► +oo. We may define a 
transfer matrix for the problem following a prescription 
proposed by Derrida, in a way to take into account the 
self-avoidance constraint exactly ||]. The connectivity 
properties of all vertical bonds of the chain arriving at a 
line ?/o from below are specified through the indication of 

1. The (unique) bond connected to the initial 
monomer of the chain (placed in y — > -co) through 
a path lying entirely below the line yo (passing 
though sites with y < yo); 

2. The pairs of bonds connected to each other through 
a path lying entirely below the line yo- 

In figure |l| the five configurations for the case m — 2 are 
depicted, and a portion of the chain placed inside the 
strip is shown with each configuration indicated. Con- 
figurations 1 and 3, as well as 4 and 5, are related to 
each other by reflection symmetry. We define column de- 
pendent activities Zi, i = 1,2, ...n a (m) according to this 
reflection symmetry, as indicated also in figure [|. The 



For a fixed connectivity configuration of a set of m + 1 
vertical bonds arriving at the line yo, the possible config- 
urations of the set of vertical bonds arriving at yo + 1 may 
be obtained, as well as the contribution to the partition 
function from the sites comprised between both sets of 
vertical bonds. This contribution will be given by 

Tla (iTl) 

i=l 

where Ni_ yo is the number of monomers with activity Zi 
in line yo and Nb iVo is the number of elementary bends 
in this line. These contributions, shown for a particular 
example in figure |^, define a line of the transfer matrix. 
For m = 2 the transfer matrix will be given by 
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FIG. 2. Contributions to the first line of the transfer ma- 
trix for the case m = 2. The connectivity configurations are 
indicated following the enumeration adopted in figure [l] 

The grand canonical partition function of the 
model, considering periodic boundary conditions in the 
y direction, will be given by 



3 = (Tr)T f 



(4) 



where N y is the total length of the strip in the y di- 
rection, so that the total number of sites is given by 
N s = N y (m + 1). The number of monomers with ac- 
tivity Zi is 



(5) 



and the total number of monomers in the chain will be 



n a (m) 

N= 

The fraction of monomers placed at column x is 

N x+ i 



p{x) 



(2-4, m / 2 )iV : 



(G) 



(7) 



where x = 0, 1, ... , n a {m) — 1 and p(x) = p(m — x). The 
Kronecker delta in the denominator contributes only for 
even values of to. 

In the thermodynamic limit N y — > 00 the partition 
function is dominated by the largest eigenvalue of the 



transfer matrix Ai, so that S ~ A^" and in this limit 



y Ai dzi 



(8) 



The first order polymerization transition in the finite 
strip will take place when the thermodynamic potential 



V> = -fc B Tln(3) 



(9) 



is equal to the one for the empty lattice ipg — 0. Thus, 
since in the thermodynamic limit we have 



i>/N s = -k B T(m + l)ln(Ai) 



(10) 



the polymerized phase will coexist with the non- 
polymerized phase for Ai = 1. Therefore, all thermody- 
namic quantities below will be calculated for Zi = z, i = 
1, . . . , n a (m) — 1 and z{n a {m)) = ojz, where the activity 
z is then fixed at the coexistence value z c so that, for a 
given value of u>b, we have Ai = 1. 

Finally, the force applied on the walls is given by 



a (<9m 



(11) 



where a is the lattice parameter and positive values cor- 
respond to attractive forces. An adimensional force per 
monomer at the coexistence may be then defined as 



/ = 



Fa _ (dm)z,u,u b _ 1 ( dz c \ 
k B TN ~ z(m) ~ z c \dm) uu ' [ ' 



Since our results correspond to integer values of m, the 
force / was estimated making the discrete approximation 

o 

f(m + 1/2, u, Lo b ) 



z c (m + 1, (j, LOb) + z c (m, ui, LOb) 
x [z c (m + 1, uj, LOb) — z c (m, ui, u>b)]- (13) 



A simple calculation may be performed in the limit 
of rigid rods u>b — 0, where bends are not allowed. 
In this limit the transfer matrix is diagonal of size 
(rn + 1) x (m + 1). For gj > 1 all monomers are on the 
walls, for to = 1 they are uniformly distributed, whereas 
for to < 1 an uniform distribution in the sites away from 
the walls is found. As expected, in this limit the force 
on the walls vanishes. In the numerical results below we 
consider u>b > 1. 



III. NUMERICAL RESULTS AND CONCLUSION 



The transfer matrices for the model were obtained ex- 
actly for strips of widths ranging from 3 to 9. After using 
the reflection symmetry, the sizes of the matrices are, re- 
spectively, equal to 6, 16, 38, 100, 256, 681, and 1805. We 
obtained the density profile at the coexistence condition 
for values of uj and u>b mostly between 1 and 3. 

For neutral walls (lu = 1) the density profile is always 
concave, with a higher density in the center. This may 
be understood since the region away from the walls is 
favored entropically. As uj becomes larger, monomers on 
the walls are energetically favored, so for sufficiently large 
values of w a convex density profile is expected. For ideal 
flexible chains (ujb = 1), at the adsorption value to = 4/3 
the density profile is flat for all sites which are not on 
the walls jl2,[L4|. Such a flat transition profile is not ob- 
served for self-avoiding chains, where convex profiles are 
separated from concave ones by an interval of values of w, 
located well below the adsorption transition value, where 
the profile is neither convex nor concave. This interval 
is quite narrow for flexible chains In figure I the 

densities are plotted as functions of to for two values of 
LOb- It is clear that as ujb is increased, favoring bends, the 
interval of values of to with a profile without well defined 
convexity grows. 
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FIG. 3. Density of monomers as function of u) for a strip 
with m = 9. (a) corresponds to uj, — 2 and (b) to u>b = 3. 
Curves are for different values of x and p(x) = p(m —x). At 
(j = 1 the highest density corresponds to the center of the 
strip (x = 4 and x = 5) and the density decreases monoton- 
ically outwards, being lowest at the wall (x — and a; = 9). 
At lj = 3 the density profile is convex in both cases, with the 
maximum density located on the walls. 

The values of uj below which the density profile is con- 
cave (wi)and those above which the profile is convex (u)2) 
are plotted in figure ^| as functions of uif, for two values 
of m. As is apparent the range with no defined convex- 
ity grows with u>t, in a nearly linear way. As a general 
rule, we found that as u)b is increased from 1, the first 
pair of densities to cross, destroying concavity, is always 
correspondent to the columns which are first and sec- 
ond neighbors to the walls. Also, the last crossing, which 
turns the profile convex, is between densities correspond- 
ing to the columns at the wall and the neighbor columns. 



In the limit io\, — > 00 we found that 0J2 ~ > 00. 
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FIG. 4. Boltzmann weights below which the density is 
concave (lower curves) and above which it is convex (upper 
curves) for m = 6 (dashed curves) and m = 9 (full curves). 

Finally, we calculated the force per monomer on the 
walls for pair of widths m, m + 1, using expression |l^. 
As already found for flexible self-avoiding chains , at- 
tractive forces appear for sufficiently large values of lj. 
Figure |5| shows results for the force as a function of uj for 
some values of u>b- The origin of attractive forces in the 
system are portions of the chain limited by monomers ad- 
sorbed on opposite walls. Thus, as expected, the curves 
/ x lu display a maximum, because the force vanishes as 
lj — > 00, since in this limit such "bridge" segments are 
absent. 




FIG. 5. Force per monomer on the walls as function of u> 
for (a) ujb = 1, (b) W& = 2, (c) cut = 3, and (d) ut — * 00. 
Results shown are for m = 8.5. 

In figure ^ the maximum force per monomer (with re- 
spect to to) is plotted as a function of cj^. It is apparent 
that the maximum of these curves is located at u>t, > 1, 
and thus we conclude that polymers for which bends are 
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somewhat favored give rise, in general, to larger attrac- 
tive forces than flexible ones. As u>b is increased, the 
maximum attractive force occurs at higher values of u> 

The force per monomer as a function of the width of 
the strip m, for uj above a threshold, shows a stable equi- 
librium point at very low separation m\ and an unstable 
equilibrium point at a larger separation 777,2 fill . The 
force is attractive in the interval m\ < x < 777,2. 
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FIG. 6. Maximum force per monomer on the walls for 
m = 6.5 (full line), m = 7.5 (dotted line), and m = 8.5 
(dashed line) as functions of Ub- Arrows indicate the maxi- 
mum force for — > oo. 

The stable equilibrium point mi(w,wt) is found in the 
range 0.5 < mi < 1.5, showing little variation as a func- 
tion of the Boltzmann factors, as long as ui is larger 
than the threshold value. The unstable equilibrium point 
7712(0;, u>b) shows a rather strong monotonic dependence 
on both variables, being an increasing function of u> and 
a decreasing function of ujb- 

In conclusion, we may summarize the behavior of semi- 
flexible chains confined inside strips observing that as the 
presence of elementary bends in the chains is favored, 
the range of values of u> for which the density profile has 
no defined convexity grows. Also, the largest attractive 
forces are found for w& > 1, that is, when the bend- 
ing of the chains is favored. It should be noticed that 
the nonmonotonic behavior observed for the tension as a 
function of the width of the strip, with two equilibrium 
points, happens in a regime where the finite size scaling 
behavior has not yet been reached |r| . 
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